One-Dimensional Impenetrable Anyons in Thermal Equilibrium. II. Determinant 
Representation for the Dynamic Correlation Functions 
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qq ' We have obtained a determinant representation for the time- and temperature-dependent field- 

f"^ , field correlation function of the impenetrable Lieb-Liniger gas of anyons through direct summation 

of the form factors. In the static case, the obtained results are shown to be equivalent to those that 
£SJ ' follow from the anyonic generalization of Lenard's formula. 
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I. INTRODUCTION AND STATEMENT OF RESULTS 

This is the second paper in a series that provides a comprehensive treatment of the properties of temperature- 
dependent correlation functions of one-dimensional (ID) impenetrable free anyons, based on the methods developed 
for impenetrable bosons The anyonic model considered in this work can be viewed as a generalization to an 
arbitrary statistics parameter k of the model of impenetrable bosons obtained from the Bose gas with repulsive 5- 
function interaction [l|, Q in the limit of infinitely large coupling constant (for other anyonic extensions of well known 
| models see 0, 5j5[l- This model, which we call the Lieb-Lini ger gas of anyons, was formulated in this form by Kundu 
. @, clarified in J 3 an d further studied in [{j [Tol. [TTI. [l2l. [l3l 1 14| . In the bosonic case, the first step in the analysis of 
the correlation functions is the derivation of the Fredholm determinant representation for these functions [15| . With 
the help of the determinant representation, a classical integrable system characterizing the correlation functions can 
be constructed as in (3,0]> leading to the short-distance and low-density expansions of the correlators. The large- 
distance asymptotics are then obtained by the inverse scattering method for the integrable system and the solution 
of the associated matrix Riemann-Hilbert problem These results will be presented in future publications. The 
purpose of this work is to derive the Fredholm-determinant representation for the temperature-dependent correlation 
— — functions in the case of anyons, and to prove the equivalence of this representation with the anyonic generalization of 
Lenard's formula [l9| . 

The results of this work can be summarized as follows. One defines free propagators 



> 

O . 

10. e(X\t,x) = e itx ~ lxX , G(t,x)= I e(X\t,x)dX, (I) 

o 
o 



and the function 



E[u\t, x) = P.V. r dX + e U\t, x)n tan ( —) , (2) 

7-oo A-/z ' V 2 / 

where P.V. denotes the Cauchy principal value. In terms of these functions, the time- and temperature-dependent 
field-field correlator of impenetrable ID anyons is: 



(^ A (x 2 ,t 2 )^ A (x 1 ,t 1 )) T = e M ^ (^G(t 12 ,x 12 ) + -^-)det(l + V T ) 



(3) 



where x a b — x a — a?6, £ a fc =ta~tb> a, & = 1, 2 , and det(l + Vt) is the Fredholm determinant of the integral operator 
with kernel 

V T (X,li) = cos 2 (™/2) exp |-^ 12 (A 2 + M 2 ) + % -x 12 {X + \J ■d(X)-d(fi) 
~E(X\t 12 ,x 12 ) - E(fx\ti, 2 ,x 12 ) a 



tt 2 {X~li) 2ir 3 



E(X\t 12 , x 12 )E(p\t 12 , x 



(4) 
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which acts on an arbitrary function /(//) as 

(V T f)(X)= I V T (\,(j,)f((j,) dfi. 



In Eq. |4|), $(\) = i)(\,T, h) is the Fermi distribution function at temperature T and chemical potential h 

®( X ' T > h ) = 1 + e (l-h)/T ■ 

Introducing integral operators Kt and A^ which act on the entire real axis and have kernels 

sinx(A — /i) 



A — fi 



and 



4(A,/i) = ^MXje Tix{X+ ^VW)- 



we obtain the static, i.e. equal-time, correlators as 

1 



(tf A (x)9 A (-x)) T = ^Tr [(1 - 1 k T )- 1 A+\ det(l - 7^T)| 7=( i +e+ — )A , 



and 



(& a (-x)*a(x)) t = ^Tr 



{l- 1 K T y 1 A^ det(l- 7 ^ T )| 7 =(i 



7=(l+e-"" t )A • 



(5) 



(6) 



(8) 



(9) 



(10) 



Here Tr [/(#, y)] = / /(a;, ar) da;, and due to the nonconservation of parity the corrrelator (^^(x)^ A (— x))t is different 
from (% A (-x)ty A (x)} T - 

The paper is organized as follows. Section |TT] introduces the Lieb-Liniger gas of anyons and presents the Bethe 
Ansatz eigenf unctions, Bethe equations, the ground state, and the thermodynamics of anyons in the impenetrable 
limit. In Section [TTT1 we compute the form factors and express the field correlator as a Fredholm determinant. Section 
IIVI presents the proof of the equivalence of Eqs. © and (fiT)|) to the anyonic version of Lenard's formula [19j. Some 
technical details of the calculations are relegated to the Appendices. 



II. THE LIEB-LINIGER GAS OF IMPENETRABLE ANYONS 



The second-quantized Hamiltonian of the Lieb-Liniger gas of ID anyons is 

fL/2 

-L/2 V 



(11) 



where c > is the coupling constant, L is the length of normalization interval, and h is the chemical potential. The 
canonical Heisenberg fields 



tf A (x, t) = e lHt ^ A (x)e- lHt , * A (x, t) = e lHt ^> A {x)e- lHt , 
obey the anyonic equal-time commutation relations 

^ A (x u t)^ A (x 2 ,t) = e- ivK < Xi - x ^<b A (x 2 ,t)^ A (x 1 ,t) + S(x 1 -x 2 ), 

* A ( Xl ,t)* A (x 2 ,t) = e <T«(-i-^ ) * A (x 2 , t )tf A (xt , t) , 



(12) 

(13) 
(14) 



^ A (x 1 ,t)^ A (x2,t) = e 



_ inKe(xi—x 2 ) 



^ A (x 2 ,t)^ A { Xl ,t) . 



(15) 



where K is the statistics parameter, which we assume to be rational (this is necessary in order for Eq. (|3T))) to hold), 
and e{x) — x/\x\, e(0) = 0. The Fock vacuum is defined as usual by 



*A(ar)|0)=0 = (0|*]i(x), (0|0) = 1. 



(16) 
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The eigenstates |Wjv) of the Hamiltonian are 

I l-L/2 rL/2 

\$n) = —= dzf- dz N xn(zi, ■ ■ ■ , z N \X lr -- ,\ N )9\{z N )~-*\(zi)\0), (17) 

VNl J-L/2 J-L/2 

where quantum-mechanical wavefunctions have the property of anyonic exchange statistics: 

X n(- ■ • ,z u z i+1 , ■■■) = e *»«(*«-*'+0 XAr (. . . , Zi+1> Zjj . . . ) . (18) 

Note that the sign in front of the statistical phase in this expression (+ittk or —ittk) depends on the choice of ordering 
of the creation operators in the definition of the eigenstates (fTT)) . The order of these operators adopted in Eq. (JTTJ) 
(leading to the phase +mitc): the particle with the first coordinate z\ created first, then z%, etc., is convenient [7J for 
the subsequent calculation of the form factors. 

In this paper, we limit our discussion to the case of infinitely strong interaction, c — * oo, which corresponds to 
impenetrable anyons. In general, the eigenfunctions xn are 0] 

XN = e + ' 2 Ej<fc£(2j Zk) J2 (_ 1)7 r e iEir =1 ^A„ (n) TTJA^ - \ n(k) - ic'e( Zj - z k )] , (19) 

sjNm j>k [{x j -x k y+^]^ 

where c' = cj cos(ttk/2), and reduce for impenetrable anyons to a simpler form: 



XN 



n^-^) e (-i)v^=^" a ^. (20) 



3>k 7tGSjv 



Here Sn is the group of permutations of N elements, and (— 1)™ is the sign of the permutation. The energy eigenvalues 

H\^ N )=E\^ N ) 
are given by the sum of effectively single-particle contributions: 



E = J2 £ ( x o), with e(X) = X 2 -h. (21) 



i=l 



The individual momenta Xj depend of the boundary conditions imposed on the wavefunctions. In contrast to 
particles of integer statistics, wavefunctions of the anyons satisfy different quasi-periodic boundary conditions in their 
different arguments, the difference resulting from the statistical phase shift 2ttk 0, S). In general, the quasi-periodic 
boundary conditions also include the external phase shift rj (we will consider rj—2'KX rational), so that the boundary 
conditions on the wavefunctions (1201) are: 



Xn(~L/2,z 2 ,--- ,zn) = e xn(L/2 1 z 2 , - ■ ■ ,z N ) , 

X n{zu-L/2,--- ,z n ) = e^-^xifizuL/^-" ,z N ), 

'■ (22) 
Xn(zi,z 2 , -L/2) = e i ^ N - 1 ^XN(zi,z 2 • • • , L/2) . 

The difference in the boundary conditions for different arguments of xn makes it possible, in general, to impose the 
condition without the statistical phase shift on any of the arguments Zj . The precise form of the Bethe equations for 
the momenta in the wavefunction (fT9|) depends on specific choice of the boundary conditions. The choice l|22[). 
in which the first coordinate z\ does not have the statistical shift in its boundary condition, gives rise to the Bethe 
equations which include the full statistical contribution ttk(N — 1) to the momentum shift of each of the anyons 
produced by the N — 1 other anyons in the system [8( : 



where r\ = rj — ttk(N — 1). Similarly to the wavefunctions, the general Bethe equations (|23p are simplified in the 
impenetrable limit c — > oo: 



e *Aji = (_l)Ar-i e iiT. (24) 
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A. Structure of the Ground State 



We assume that the ground state of the gas contains N anyons, and take, for convenience, N to be even, although 
this does not affect our final results. We denote the momenta of the particles in the ground state as fij, where 
j = 1,- ■■ , N, and introduce the notation {[...]} such that 

{[ar]}=7, if x = 2tt x integer + 2n^ , 7 6 (—1,1). (25) 

The Bethe equations (|24p give then the momenta fXj : 

2tt / N + 1 \ 2tt<5 
^ = Ty^~ ) + ~L~ ■ ^ = 1 '"-' iV 0' ( 26 ) 

where 5 = {[rj]}. In the thermodynamic limit L — > 00, N ~ * 00, N/L = D, momenta of the particles fill densely the 
Fermi sea [—q, q], where q — yh is the Fermi momentum and the gas density is D = q/n. 



B. Thermodynamics 

The thermodynamics of the Lieb-Liniger anyonic gas was considered in [l(| [ll[ . Similarly to the structure of the 
ground state, all local thermodynamic characteristics in the case of impenetrable anyons are equivalent to those of 
the free fermions. At non- vanishing temperature T, the quasiparticle distribution is given by the Fermi weight ((6|), 
and the density and energy are 

1 roc 1 rco 

D=— d{X,h,T)d\, E = — \ 2 $(\,h,T) dX. (27) 
2^ /-oo 2tt J_ qo 

The density increases monotonically as a function of the chemical potential h. At T = 0,we have D = for h < 0, 
and 0<-D<ooifO<ft.< 00. At non-vanishing temperature, the density is zero for h = —00 and monotonically 
increases with h for — 00 < h < 00. 



III. TIME DEPENDENT FIELD-FIELD CORRELATOR 

In our previous paper 19], we have derived the anyonic generalization of the Lenard formula, which for impenetrable 
free anyons, is an expansion of the anyonic reduced density matrices in terms of the reduced density matrices of free 
fermions. In the simplest case, the correlator 

(zi|p?|z 2 ) - (^ a (x2)9 a ( Xi ))t (28) 

is the first Fredholm minor of an integral operator, whose kernel is the Fourier transform of the Fermi weight ©. In 
this section, we obtain the time dependent generalization of this result. Our approach will be based on the following 
considerations. We start with the zero temperature field correlator 

(^ A {x2,h)^ A {x 1 ,t 1 )) N = — — , (29) 

(V(/ti,- ■ • ,fl N )\W(fli,- ■ ■ ,flN)) 

where the wavefunctions are taken to be normalized as 

(*(mi,--- ,Vn)) = L n , (30) 

and hi, ■ ■ ■ , hn are the momenta in the ground state (f2"o| . Using the resolution of identity for the Hilbert space of 
N + 1 particles 

* = V- |^(Ai,--- ,Aat + i))(^(Ai,--- ,A w+1 )| 

fj (*(A 1 ,---,A W+1 )|*(A 1 ,--- ,X N+ i)) ' [ > 

all {A} N+ i 

where, according to ([30]) 

(*(Ai, ■ • • , X N+ i)\*(Xi, ■■■ , Ajv+i)) = L N+1 , 
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and the sum is over all possible solutions of the Bethe equations with N + 1 particles, we have 

(*A(s3,fc)*Ji(zi,ti)>JV = I ^ TT E (*Ar({M»l^(a:2,t a )|* W+ i({A}))<* w+1 ({A})|*] 1 ( a! i ) t 1 )|* Ar ({/i})) . 

all{A} JV+ i 

(32) 

Dehning the form factors 

F N+ltN (x,t) = (^ +1 ({A»|*1 1 (a:,t)|* JV ({/i})) , i^ +liJV 0M) = (* JV ({ y u»|^ j4 (c C , *)l*w+i({A») , (33) 
we can rewrite Eq. (|32[) as 

(^(3:2,^2)^(3:1, ti))jv = J2W+1 F N+l,N( x 2,t 2 )F N+1 . N (x 1 ,t 1 ) . (34) 

all{x} N+1 

Equation (|34[) means that in order to find the dynamic field correlator, we need to compute the form factors and 
sum over all of them. After the summation, one can take the thermodynamic limit. In general, such a summation 
of form factors is extremely difficult. The main simplification which makes it possible to perform this summation 
in the model of anyons we consider here, is the fact that, similarly to the problem of impenetrable bosons [H, [Til ], 
the local thermodynamic properties of particles with ^-function interaction are identical with those of free fermions 
regardless of the actual exchange statistics. Finally, the finite-temperature correlator can be obtained from the zero- 
temperature result using the standard argument developed for the Bose gas (see, e.g., Appendix XIII. 1 of [l|), which 
is also applicable in the case of anyons. 



A. Form Factors 



As a first step in carrying out the program outlined above, we compute the form factors. In the definition (f3"3")) of 
the form factors, the eigenstates |^jv({/x})), |\E r jv+i({A})) have, respectively, N and N + 1 particles. Although the set 
{^} represents in (f3"3"| momenta in the ground state of N particles, our calculation below is valid also when \^ N ({^i})) 
is not the ground state. As before, we assume for convenience that N is even. We denote by {/ij} the momenta of 
the anyons in the A''- particle eigenstate, and by {\j} the momenta in the N + 1 eigenstate. 

Using the definition (fT"7|) for the eigenstates with N and N + 1 anyons 

I*aKM)> = "4f / d N z X n(zi, ■ • • , zn\{h})& a {zn) ■ ■ ■ &a(*i)\0) , 



(**+i({A})| - -/== J dN+1 y (o|*a(i/i) • ■ ■ Mwv+ikWw. • • • . y^+il{A}) 

one can write the form factor as 

F N+1:N (x,0) = -y= = == J d N+1 y d N z x*N+i(yir ' ' ■ Vn+i\{>})xn(zi,- • • , z N \{fi}) ■ 

(0\Va(vi) ■ ■ ■ *A(yN+i)* A (x)& A (z N ) ■ ■ ■ & A (zi)\0) . (35) 

A direct application of the anyonic commutation relation (j 13|) and Eq. (|16[) described in more details in Appendix [XJ 
reduces this expression to 

0) = {9tf +1 \^ A (x)\9tf) = VATT / d N z x*n+i(*i, ■", zn,x\{X})xn 



An important feature of Eq. (f3T>)) is that the order of the creation operators chosen in Eq. (| 1 7[) makes the "free" 
coordinate x in (|36p the last argument of the wavefunction xn+i- This ensures that both wavefunctions, xn and 
Xn+i, have the same phase shifts (|22[) at the boundary of the normalization interval in all other variables Zj. Since 
these phase shifts are canceled in Eq. , the expression under the integrals over zj is periodic in each of the variable 
Q. This feature is the necessary consistency condition for the Hilbert spaces of anyon wavefunctions with different 
numbers of particles, and is important in what follows for the appropriate calculation of the form factors (|36|) . 
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The sets of momenta {/Zj} and {Xj} in the wavefunctions xn and xn+i in (]36|) are determined by the Bethe 
equations (fM]) as 

Mj = ?E f mj + V\ + 5 = {[t 1 -ttk(N-1))}, j = l,--. ,N, mj eZ, (37) 

2tt 2it5' 

X i = T nj + ~' 6 ' = {[V-^N]}, j = !,■■■ ,N+1, nj £Z. (38) 



These equations show that 



which means that Aj and [ik never coincide except in the trivial case k = 1, when we have a gas of non- interacting 
fermions. In all other situations, Xj and /i/. are different. This difference between them comes from the phase shift due 
to the hard-core condition on the added particle described by the factor 1/2 in (|3TH) , and the extra anyonic statistical 
phase added to the anyon system together with the particle [7j . This difference between Xj and /ifc plays an important 
role in the following calculations. Using the identity 

e +t ^ e{x - v) e(y - x) = cos (ttk/2) e(y - x) - i sin (ttk/2) , (40) 

we can rewrite the anyonic wavefunction (f2"0"]) as 

FT [cos (ttk/2) e(z.j — zu) — isin (7r«;/2)l . 

X ^l,-, Zjy |W)^ llj>U vW! £(-l)* e »52=X*»*«0. (41) 



Using this expression for both of the wavefunctions in (|36p we obtain 
Fjv+^ar.O) = -1 ]T £ (_i)^ e -^ ( „ +1) 



iV! 



X 



ri/2 JV_ 
'-V2 „=l 

Integration by parts in this equation produces the boundary terms in the following form 



/■i/2 N 

/ FT dz„ [cos (ttk/2) e(x - z„) + i sin (ttk/2)] e _i S»=i *«( a «c»)-^(«)) . (42) 



— — (cos (ttk/2) e(x ~ z n ) + i sin (ttk/2)) 

-H'W(ri) - M<r(n)J 



=L/2 



-L/2 



2 e 2 



*7(A„(„)-At ff(n )) 



*(^7r(n) ~ Mff(n)) 

All these terms vanish due to Eq. (|39| . Then, using the relation 

de(x - z n ) 



_|_ g+"TK; e iL(A jr ( n) -/V(»))j 



we obtain the following expression for the form factors 

N 



= -26(x ~ z n ) , 



. [2*0)8(^/6/2)]^ 

Fjv+i,/v(a;,0) = - v exp ^ /.r 



JV+l 

53 ^ - 53 a j 



JV 



E Ec-^Ih — - — 



(43) 



(44) 



(45) 



This expression differs from the corresponding result for impenetrable bosons P,[lB] by the spectrum of the momenta 
which now include the statistical shift, and by the overall [cos(ttk/2)] n factor. For k = 0, both differences disappear, 
and Eq. (|43]l reproduces, as should be, the case of the Bose gas. We transform this equation following the corresponding 
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steps for bosons [l|, [H[ . One can see directly that the sums over permutations in (|45f can be written in the form of 
a determinant: 



m E Enrn^ 



=1 ^r(i) - /Mi) 



= ( 1 + 7?- IdetvfM. 



9a 



a=0 



with 



1 



Xj — fik Xn+i — M/c 



, 3,k = l,--- ,N, 



(46) 



(47) 



reducing Eq. (|45j) to 

^V+i,jv(^) 0) = (2i cos(ttk/2)) n exp t ix 



N 



N+l 

E^ _ E A i 



detjy (M; fe ; 



(48) 



a=0 



The determinant part of this equation can also be written as 

d 



1 



dc 



det N {Mf k ) 



ttES 



E HO* 11^ — 

j=i X «U) - h 



(49) 



as one can see directly from the L.H.S. of (|46|) by noticing that due to the permutations tt of Xj, all permutations of 
fij give identical contributions to the sum over tt € Sn+i- 

Alternatively, one can introduce a fictitious momentum /ijv+ij and obtain the following representation [20| of the 
form factor in terms of this momentum: 



Fn+i,n(x,0) = (2i cos(7tk/2)) exp < ix 



N 



N+l 

E^' - E a j 



lim 

fwv+i— >°° 



-fJ>N+i detAr +1 



1 



(50) 



where detjv+i(%'fc) is the determinant of the (N + 1) x (N + 1) matrix with elements Oj&. We will not be using this 
representation explicitly below. 

The time-dependent form factors can be obtained from the timeless form (|48[) using the following simple relations: 



iHt 



* N m)) = e- it ^=^-V\* N ({ri)) 



and 



(* N ({X})\e im = e lt ^^- h \^ N ({X})\ 



(51) 



(52) 



Combining the exponential factors in these expressions with those in Eq. (|48[) . we arrive at the final result for the 
time-dependent form factor: 



/N+l 



N 



F N+1 , N (x,t) = (2icos(n K /2)) N e- M ( J] e(X t \t,x) j [ J] e*{^\t, x) ] (l 



9a 



det^v (M? fc ) 



a=0 



where we have introduced the function 



e(A|t, x) = e 



(53) 



(54) 



e*(X\t, ir) is its complex conjugate, and Mjj. is defined in (j47|) . The form factor of the annihilation operator ^>a(x, t) 
is obtained through complex conjugation 



(*Jv(M)l*A(a:,*)|*Ar+i({A})) = Jft +1 , w (s,t) ■ 



(55) 
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B. Summation of the Form Factors 



Using Eqs. (|53|) and |55|) , we write the field correlator (f34| as a sum over intermediate momenta {A}: 



all {a}jv+i 
d 



(2co£(W2)r 

L2N+1 



'N+l 



N 



J\ e*(Xi\t 2 i,X2i) J JJ e(fij\t2i, X21) 



X 1 + 



<9f 



det^v {Mf k ) 



a=0 



. i=l 

H 0/3 



det 



K) 



(56) 



0=0 



with the notations a; a fc = x a — Xb, t a {, = t a — %, a, b = 1, 2. The matrix M jfe here is the same as (|47p with a replaced 

by /3. As was mentioned above, modulo the [cos(7TK/2)] 2Ar factors and the spectrum of momenta, Eq. (|56|) is identical 
with the expression for the bosonic field correlators [J [HI]. This means that the summation process over {A} is very 
similar, and we just sketch the derivation here. Since we sum over all momenta {A}, individual momenta Xj are 
equivalent up to permutation. This means that one of the two permutations of {Xj} involved in the definition of the 
two determinants in (|56[) produces coinciding terms, so that under the sum over {Xj}, one can replace one of the 
determinants, e.g., the second one, with 



N 



(57) 



obtaining 

{^ A (x2,t 2 )^ A (x 1 ,ti)} N 

= e lht ^ \\[e( N \t 21l x 21 ) 



1 ( 2cos(7tk/2) 
L 



L 



2N 



x E 

all {A}jv-fi 



d 



e*(Ajy+i \t 2 i,x 2 x) + — ) det at 



(N + l)\ 

e*(Xj\t 21 ,x 2 i) 
(Xj - fi k )(Xj - Hj) 



e* (A,- |t2i,a:2i) e*(X N+ i\t 2 i,x 21 ) 

t 

(Xj-Hj) (X N+ i-(ij) 



(.58) 



The summation over the momenta {Xj} can be done then independently over each A^ inside the determinant. Also, 
we transfer the factors e(fj,j\t 2 i,x 2 i) in (|58p into the determinant splitting them between the rows and columns, and 
use the formula 



1 



1 



1 



1 



(Xj - ilk) (Xj - Hj) \Xj - (j,j Xj - fi k J fij - ii k 



1 



(59) 



This gives the correlator as 



(*A(x2,t 2 )* A (xi,h)) N = e lht ^ (^-G L (t 12 ,x 12 ) + — 

\2tt oa 



x det 



N 



x / 2(1 -d jk ) 
\irL(^j,j - fi k 

where we have defined the functions 



Sj k E L ([j, k \ti 2 ,xi 2 )e([ij\t 21 ,x 2 i) + e(^ J \t 2 ,x 2 )e*(n k \ti,x 1 )cos (ttk/2) 

a 



\ ~~T7 i ^ ( E L (Mi I h2,xi 2 ) - E L (n k ,ti 2 ,xi 2 )) - —^E L (iJ,j\t 12 ,x 12 )E L (n k \ti 2 ,xi 2 

\irL(fj,j - fi k ) L-k z 



a=0 



±-G L (t,x) = ±J2< X \ t > x )> 



^w>*) = zE 



e(X\t, x) 



L X - [i k 



(60) 

(61) 
(62) 



E L (Hk\t,x) 



4cos(ttk/2) 2 
L2 



E 



e(X\t, x) 
(A-Mfc) 2 ' 



(63) 



and A = ^£(Z + 5') - see (j3HJ). Formula §^ is the final expression for the field correlator in the ground state of N 
anyons on a finite interval with quasi-periodic boundary conditions. 



9 



C. Thermodynamic Limit 



In order to obtain the correlator in the thermodynamic limit, we need to compute the large- L limit of the functions 
([6T]). and ((Ml)- This is done in Appendix [B] with the results 



/OO 
e(X\t,x) dX, 
-oo 

E(ji k \t,x) = lim E L (n k \t,x)=P.V. [ 

L — oo J_ 



,, e(X\t,x) . . . / 7TK 

dX — 1- e(/Lt fe i,x)7rtan — 

A - v 2 



~. , . 1 , ~ , . . /,N 2cos 2 (7tk/2) 9 

E{fj, k \t,x)= hm E L {fj, k \t,x) = e{n k \t,x) H — E{fj, k \t,x). 

L^oo irL OjjL k 



(64) 
(65) 
(66) 



In the thermodynamic limit L, N — > oo with D = N/L constant, the anyon momenta fill densely the Fermi interval 
[— q, q], where q = \fh and D = q/ir. In this case, the determinant in the correlator (pO)) can be understood as 
the Fredholm determinant of an integral operator. Indeed, for an arbitrary integral operator V, whose action on a 
function /(A) is defined by 

(67) 



(Vf)(X) = J V(\,fi)f(ji) dpi, 
the associated Fredholm determinant is (see, e.g., |22j ) 



det(l + V) = lim 



l + £V(Ai,Ai) £V(Ai,A 2 ) 
£V(A 2 ,Ai) i + s;v(x 2 ,x 2 ) 

£V(X n ,Xi) £V(X n ,X 2 ) 



£V(Xi,X n ) 

av(x 2 ,x n ) 
i + iv{x n ,x n ) 



(68) 



where £ = (b — a)/n, X p — A p _i = £ and Ao = a, A„ = b. One can see directly that, in the thermodynamic limit, the 
determinant part of Eq. (|60|) has the same structure with N momenta fij separated by £ = 2ir/L filling the Fermi 
interval [—q, q]. This means that the correlator can be expressed as 



a{x 2 M)^\{xiM)) = 



where Vq acts on an arbitrary function /(A) as 



^G(t 12 ,x 12 ) + -^),hu/{ +]],< 



a=0 



(Vof)(X)= I F (A, m )/(m) dn, 



and 

Vb(A,/x) = cos 2 (7TK/2)e(A|i 2 ,a;2)e*(^|ii,a;i) 
Performing the unitary transformation 

V (X,fi) = exp <{ -i 
with the property 

det(l + Vo) = det(l + %) , 
we transform the kernel Vq(X, fi) (JTTJ) into the symmetric form: 



E(X\t 12 ,x 12 ) - E(n\t 12 ,x 12 ) a . 
tt 2 (A - fi) ^E(X\t 12 , x 12 )E(p\t 12 , x 12 ) 



. (h + t 2 ) 2 2 

(A — /x ) + i 



{xi + x 2 ) 



(A -a*) Vb(A,M), 



(69) 



(70) 



(71) 



(72) 



(73) 



V (A,/i) = cos^ttk/2) exp <^ --*i 2 (A 2 + //) + ^x 12 (A + /*) 



7r 2 (A-/i) ^£(A|tl 2 ,Xl 2 )£(Ml*12,*l 2 ) 



(74) 
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Two observations are in order. First, one can check that the second term in (|66|) is obtained from the first term in the 
square bracket of (f?4")) in the limit A — ► fi. Second, in the limit k — ► 0, Eq. (|74[) reproduces the known result [3, [l5| 
for impenetrable bosons. 

In the static case (ti = t 2 ), which is discussed in the next Section, the kernel (f?4")) can be simplified further. One 
needs to distinguish two cases. 

• x\ > x 2 ■ In this case, 

E(X\0,x 12 ) = -ine- iXl2X [l + itan(Trn/2)] , (75) 

and the kernel (1741) becomes 

(76) 



• xi < x 2 . In this case, 

E(X\0,x 12 ) = i7re- ixi2A [l-itan(7rK/2)], (77) 

and 

T/-a \ (l + e— ) f S in(z 12 (A-A0/2) ^ q _„ K r .x 12 o 

y (A,„) = ^ — j+-e expj-z — (A + /x)j. (78) 

We now extend the discussion to the situation of non-vanishing temperature T. The temperature-dependent field 
correlator is defined as 



Tr 



{e- H / T ^ A (x2,t 2 )^ A (x 1 ,t 1 ) 



(* A (x 2 ,t 2 )* A (x 1 ,t 2 )) T =^ Tr e _ H/T '-■ (79) 

According to the well-known argument developed for the Bose gas [l[ , this correlator can be found as the mean value 
over any one of the "typical" eigenfunctions fix of the Hamiltonian which characterizes the given state of thermal 
equilibrium: 

{n T \^ A (x 2 ,t 2 )^ A (x 1 ,t 2 )\n T ) 

MM ■ (80) 

This argument depends only on the general saddle-point approximation in the description of the state of equilibrium, 
and also holds in the case of anyons. The further computation of the field correlator based on Eq. ([80)) is similar to 
the zero-temperature case, the main difference being the change of the measure of integration: 



dX I dXti(X,T,h) with d(X,T 7 h) = * , WT , . (81) 

) 1 + e^ A n 'i 1 



The final result for the temperature-dependent correlator is then 



{■$ A (x 2 ,t 2 )V\(x 1 ,t 1 )) T = e l,lt21 ( ^-G{t 12 ,x 12 ) + 4z )det{l + V T ) 



(82) 



a=0 



where the kernel of the integral operator Vr is 



v T (x,») = ^MXjv Q (x,f,)yMri, 

2 (7TK/2)exp|-^ K 
' E(X\t 12 ,xi 2 ) - E(p\ti 2 ,x 12 ) o_ 

7T 2 (A — jU) 27T 3 

and the operator acts on an arbitrary function f(/i) as 



cos 2 (™/2) exp <j --t 12 (X 2 + fi 2 ) + -x 12 (X + fi)j y/0(\)0(jj,) 

E(X\ti 2 , xi 2 )E(p\ti 2 , xi 2 ) 



(83) 



(V T f)(X)= / V T (X,fx)f( t x)dfi. (84) 
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IV. EQUIVALENCE WITH LENARD FORMULA 

In the earlier paper [l9l |. we obtained the anyonic generalization of the Lenard formula for the equal-time field 
correlator or, equivalently, reduced density matrices of anyons. In the case of the first reduced density matrix, the 
anyonic Lenard formula reads 



(x\fi\x') 



where the kernel of the integral operators 8 T is 



det 1 - 70= 



7=(l+e±""')/7r 



M£ - V) = 5 



,»U-i7)A 



l + e (X>-h)/T ' 



and their action on an arbitrary function is defined as 



(<?t/)(0 = 



H - V)f(v) drj. 



(85) 



(86) 



(87) 



i± 



In these expressions, the plus sign refers to the situation when x' > x and 7+ = [x, x'], and the minus sign - to 
the situation when x' < x and J_ = [a:', a:]. The resolvent kernels associated with the kernel 9t(x,u) acting on the 
intervals I± are denoted by (£, 77) and satisfy the equations: 



(1 + e ±llXK ) 



Ji± 



One can rewrite Eq. (|85j) in terms of the resolvent kernel qt and the field correlator as [If 

{& a (x')*a(x))t,± = ~Qt(x', x) det (l - 7 0±) 



7=(l+e± i ™)/ 7 r 



(89) 



where again, the plus sign refers to the case x' > x and the minus sign - to x < x' . Next, we show that Eq. (|89p is 
reproduced by the results obtained in the previous section when they are specialized to the equal-time correlators. 
We treat the two cases, x' > x and x' < x, separately. 

A. The static correlator {9 A (-x)^f A (x)) T 
Equations (fB"4")) and ([M]) show that in the static case 



—G(0,x)=S(x). 
Zn 

Using this relation and Eqs. (|T6"|) and (f8"3")l . we see that the equal-time field correlator can be written as 

(9 A (-x)tf A (x)) T = \ 5(2x) + t^J det f 1 - { — >-K T + a—A+ 

where Kt and A^, are the integral operators acting on the real axis and defined by kernels 

sinx(A — /i) 



K T (\,ri = y/0(Xj 



X — (i 



and 



At zero temperature, both operators act on the interval [—q, q] and their kernels are 

X — ji 



(90) 



(91) 



(92) 



(93) 



(94) 
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The commutation relation (|13p shows that 

(* A (-x)tf A (x)) T = e™ K (* f A (x)* A (-x)) T + 6(2x) . (95) 
This means that in order to prove the equivalence with Lenard formula, we have to show that 

= a t *»(tf A (x)9 A (-x)) T) (96) 



G + (k,x,T) = ^-dct — L Kt + 01 ^ A + 
dot \ 7r 2tt 



a=0 



where {^ A (x)^! A (—x))x is given by (|89|) . For a general integral operator with kernel V, one of the useful expressions 
for the Fredholm determinant is 



In det(l - jV) = - V — Tr V 
^-^ n 



n=l 

Making use of this formula, we obtain 



G+(k,x,T) = — Tr [(l- 7 K T )- 1 i+jdet(l-7^T)| 7 =( 1+e *-)/, r . (97) 



Denoting as /_ (A) the solution of the integral equation 



f i J- P i7rK \ r°° , 

/+( A ) _ L_ i / K T (X,n)ft(j*) dfx = v^(A)e-" A , (98) 



we can rewrite ([§7]) as 



G+(k,x,T) = —J e-™ x f^(\)^d(Xjd\ det(l- 7 ^r)| 7 =(i +e <-«)/ ff . (99) 
We will show now that 

det(l - jKt) = det (l - 7 0+) , (100) 

where the operator (9<r is described by Eqs. (|56")l and ([57]) . and 7 = (l + e l7rK )/7r. Direct and inverse Fourier transforms 
of a function g can be defined to include as integration measure yj $(X): 

-1 /"OO />0O 

g(A) = / d£ e^ ff (0 , = / dX ^{X)e-^g(X) . (101) 

With this definition, taking the Fourier transform of the integral equation 

g(0 - 7 T M£ - f MO # - G(0 , (102) 

*/ —x 

we obtain 

/•OO 

g(A) -7/ K T (\-p)g{p) dn = G(X). (103) 



Coincidence of the two equations implies the equality (| 1 00|) of the determinants. 
The final step in proving the equivalence of Eqs. (|89|) and (|9T|) is to show that 

= i |°° e-" x fl{X)y/d{X)dX. (104) 

The Fourier transform of the equation defining the resolvent kernel qt 

(1 4- e inK ) f x 

<&(Z,-x)-- -/ 0T(t-?)(&(t',-x)dS' = 6T(Z + x), (105) 
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gives 

~ Q +(\-x)- y — '- K T (X-^g+(p,-x)dfi=-e- lxX ^{X). (106) 

7T J-eo 2 

Comparison of this equation with the definition of /_ (A) shows that 

~Q+{\,-x) = \f+{\). (107) 

Taking the inverse Fourier transform of (|107|) proves (|104[) . Thus, we have shown that for x' > x, the Lenard formula 
([55)) is equivalent with the result ([51]) for the static field correlator that follows from the direct summation of the form 
factors. 



B. The static correlator (* a(x)^ a (-x))t 

In this case, the proof of the equivalence of the two approaches is very similar to what was just discussed for x' > x. 
Equations (|78j) and (f83|) show that the static field correlator is 



(108) 



a=0 



(* A (x)tf A (-x)) T = U(2x) + £J det M - [ — - >-K T + ^ T 

where Kt is given by (|92|) and 

A T (A, fi) = ^/J(Xje^ x +^ yfiffi . (109) 
From the commutation relation (fT3"|) we see that 

(9 A (x)tf A (-x)) T = e^^-x^^T + S(2x) , (110) 



so we have to show that 



d ( (1 + e~™ K ) ~ 
G-(k,x,T) = — det 1- ±— >-K T + a AZ 

OOi V 7T 2tt 



= e-™(* A (-x)* A (x)) T . (Ill) 

a=0 



where (^ A (— x^a^^t is given by Eq. (f8"9"|) . Similarly to the discussion in the previous section, we can rewrite G 
as 

-ittk poo 



G~{n,x,T) = e — — / e +ixX f+{X)^d(X)d\ det(l - 7 ^r)| 7 =(i +e -^)/ w , (112) 
where /7 (A) is the solution of the integral equation 

f 1 J. p 47TK\ POO 

/+ (A) — ^— — '- K T (\^)f-(n)dv=VW)e +rxX - (H3) 



The equality of the Fredholm determinants of the operators Kt and 9t was shown in the previous Section, so it 
remains to prove that 



i r°° , 

(-x,x) = - e + " A ^(A)v^(A)rfA. (114) 



g T v -,-v 2 



Again, taking the Fourier transform of 



we obtain 



which shows that 



Cl 4- e~ i7rK ) 

-/ M^-Oer(^.«)de' = Mf-a;). (115) 



(116) 



g-(A,a:)-^ '- / X T (A-^)g T (^,a ; )d M = -e+" A v ^(Ay, 



~g T (\,x) = ±f+(\). (117) 



The inverse Fourier transform of (|117|) gives the correct result (|114p . 



14 



V. CONCLUSIONS 



In summary, we have obtained the time- and temperature-dependent correlation functions of fields for impenetrable 
ID anyons as Fredholm determinants. The Fourier transform of the corresponding integral equations proves the 
equivalence of our approach with the anyonic Lenard formula derived previously (Eq. 57 of [l!|) for the one-particle 
reduced density matrix of anyons. The same technique can be used to obtain the multi-point correlation functions 
from the Lenard formula for n-particle reduced density matrices (Eq. 56 of p^|). The next step in the exact calculation 
of the anyonic correlation functions is to use the determinant representation derived in this work to obtain a classical 
integrable system of nonlinear differential equations characterizing these functions. These equations should make it 
possible to construct the short-distance and low-density expansions for the correlators. This will be addressed in a 
future publication. 
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APPENDIX A: ANYONIC FORM FACTORS 



In this appendix, we prove Eq. ([36|) . Consider first the simple example of the form factor F 3i2 : 
1 



*3,2(a:) 



d 3 y d 2 Z X3(yi,2/2,y3)X2(^l,22)(0|*A(yi)* J 4(2/2)*A(j/3)^! 1 (a;)^(z2)*]4(^l)l > 



If one defines 

A = (0\^ A (y 1 )^ A (y 2 )^ A (y 3 )¥ A (x)^\(z 2 )^ A (z 1 )\0) , 
then successive applications of the commutation relation (fT3")) followed by the Eq. (1 16[) give 

A = (0\y A ( yi )* A (y 2 ) [y A (x)y A (y 3 )e-^<y^+6(y 3 ~x)] ^ A (z 2 )M> A ( Zl )\0) 

= (0\* A ( yi )V A (y 2 )* A (x) [f A (z 2 )^A(y3)e~ t7rKe{y3 - Z2) + 6(y 3 - z 2 )] ^ A (z 1 )\0)e- i " te ^- <e ^ 

+(0\^ A (y 1 )^ A (y 2 )^ A (z 2 )^ A (z 1 )\0)S(y 3 - x) 
= (0\* A (yi)* A (y 2 )* A (x)* A (z 2 )\0)5(y 3 - Zl ) e -*»«M»»-*0+"(v»-*>] 

V v ' 

(a) 

+ A { Vl )^ A {y 2 )^ A {x)^ A {z 1 )\Q)5{y 3 - z 2 )e~^ ^-*) 

V V ' 

(b) 

+ {Q\^ A {y 1 )^ A (y 2 )^ A {z 2 )^ A {z 1 )\Q)5(y 3 - x) . 



(c) 

Performing similar transformations, we obtain 

a = %i - x)S(y 2 - z 2 )S(y 3 - « 1 ) e -*' r ' t Wva-*)+=(»8-« a )+e(i(8-x)] 
+S( Vl - z 2 )5(y 2 - x)S(y 3 - Zl ) e -""*(!/3-* 2 )+«(v3-*)] , 

b = 8{y x - x)5{y 2 - Zl )5(y 3 - ^"^(^H^-*)] 
+S( yi - Zl )5(y 2 - x)S(y 3 - z 2 ) e — ( 

c = S( yi ~ z 2 )S(y 2 - Zl )5(y 3 - x ) e ^^'^ + 5{ Vl - Zl )S(y 2 - z 2 )8{y 3 - z 3 ) . 
Substituting A = a + b + c into (|Aip , we have for the form factor 



*3,2(a0 



l -j= f d 2 Z {x3(^^2^l)X2(^1^2)^^ [£( ^^ )+e(Zl " Z2)+e(Zl -^ )1 



2V3 

+X* 3 (z 2 , x, Zl ) X 2(zi, Z2 )e-^<zi-*2)+<^-*)} + x *( Xj zu z 2 ) X2 ( Zl , « 2 )e-*T«[e(*i-)+<» a -*)] 
+Xl(zi,x, z 2 )X2(zi,z 2 )e- im ^-^ + xl(z2,zux)x2(zuz 2 )e~ 1 ^-^ 

+X3(2l; Z 2, X)X2(Z1, Z 2 )} . 



(Al) 
(A2) 



(A3) 



(A4) 

(A5) 
(A6) 



(A7) 
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Using the anyonic property (|18p of the wavefunctions, and its complex conjugate: 

X *(... ,z h z i+ i r --)=e- i ™^-*^ X *{--- (A8) 
we reduce Eq. IA7I to the final expression for the form factor 

Fz, 2 {x) = VS J d 2 zxl(zi,Z2,x)x2(zi,z 2 ). (A9) 

The calculations leading to Eq. (|A9|) can be generalized to arbitrary N: 

F n +i,n(x) = {^ N+ i\^ A (x)\y N ) = VN+1 J d N z x*n+i{zi,--- ,z n ,x)xn(zi,--- ,z n ). (A10) 

This result follows from Eq. (|35|) by noticing that the statistical phase factors in the commutation relations (fl3]) 
(fT5")) of the field operators are compensated by the exchange property (fT5|) of the wavefunctions. This means that 
the pairing of the ^^(x) operator with any of the 'J Aijjj) operators produces N + 1 identical terms in which the 
coordinate x is made the last coordinate of the wavefunction xn+i- After that, the integrals over z's and remaining 
y's can be limited to the ordered regions z\ > z 2 > ... > Zn and yi > y 2 > ■■■ > UN giving directly (|A10[) . 



APPENDIX B: THERMODYNAMIC LIMIT OF SINGULAR SUMS 

In this appendix, we study the behavior of the functions defined by Eqs. (|6Tj) . (162|) . and (|63|) in the thermodynamic 
limit of large length L of normalization interval. We start with In this case, the function summed over the 

momenta A is sufficiently smooth, so that the anyonic shift 2ttS' / L of the momenta becomes negligible when L — > oo, 
and one can pass directly from the sum to the integral over A: 

o />oo 

G(t, x) = lim G L {t,x) = — V e(X j \t,x) = e{X\t,x)dX. (Bl) 

L— »oo L z — ' / 

The regularization t — » t + iO for e(X\t, x) — exp(ziA 2 — ixX) is implied in these expressions. 

Next, we turn to Eq. (|62p . In this case, the function under the sum is no longer smooth in the thermodynamic 
limit. We transform it by separating the singular part that can be summed explicitly: 

i i . ,.„/■,■> 27r e(\j\t,x) 

E(ii k \t,x)= lim E L ( P k\t,x) = — V ; 3 1 ' 

L^oo L L — ' A, — lit 

2tt v-* e(Xi\t,x) — e(u,k\t,x) . , , / k + 1\ 1 ,_ . 



a 3 . e¥(z+*') 



n— — oo 



In the last line here we have used Eq. The first term in (|B2[) is now a smooth function, so as before, we can 

directly replace the sum with the integral, since the anyonic shift of the momenta does not affect the value of the 
integral. The integral can then be transformed as follows: 



dX <Mt,x)~e(, k \t,x) = Ry _ r°° dX ej^) Ry< r dX 



A - Mfc J-oo A - ^fc 7-oo A - //ft 

= P.V.fdA«^. (B3) 

./-oo X — (Ik 

Under the natural interpretation of the sum in the second term in (|B2[) . it can be simplified using formula 1.421.(3) 
of [Hj], 7rcot(7ra;) = (l/x) + 2xJ2™ =1 {x 2 ~ n 2 )- 1 : 

n— — oo v 7 

Collecting the two terms we finally get 

Efr k \t, x) = P. V. / dX e . (A| * ,x) + e(/i fc |t, x)tt tan . (B5) 

J-oo A - Mfe V ^ 7 
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The function denned by Eq. (|63[) is more singular than E(nk\t, x) f 1|B2|) . To transform it, we use the same strategy 
of separating the most divergent terms that can be summed explicitly: 

~ , 4 o, , . v—v e(\j\t,x) 

E (n k t,x = lim E L (ji k \t,x = -^cos 2 ttk/2 V y Jl ' , 

L^oo L z * — ' A,- — Ilk) 

^ cos2(W / E e feM)-eWM) +e(wM) ^g _. y (B6) 

/(Mfe)= ^ e (A 3 |t,x)- e ( Mfc |t, 3 :) > (B7) 

In,- V "^ l '- r '- M/ 'f l '- r, = -^(p.V. / dA^^j+^^.tanf-) . (B9) 

For the second term in the R.H.S. of (|B6]) we use the formula 1.422.(4) of [H[ tt 2 / sin 2 (7nz;) = J2n=-oo( n ~ x )~ 2 to 
get 

«+l^- 2 



Defining 



one has 



Taking the limit L — > oo and using (|B3[) and (|B4[) in this equation we obtain 

L ^ e(A i |t, x) - e(^ fc |i,x) _ 5 / f 00 e(A|i,a;)\ 9e(/i fe |t,a;) 



71 — — CO 



("- — ) =sp)' (B10) 



Collecting all the terms we have the final result 

s, , , v , 2cos 2 (7r K /2)de(^M) /tt/cn 2 cos 2 (7m/2) 9 / [°° e(\\t,x) \ 

E(nk \t, x) = e{n k t, x) H = Trtan I — M — P.V. / dX - . (Bll) 
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